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Abstract 

Constraints on the two Higgs doublet model are presented, assuming a top 
mass of 174 ± 17 GeV. We concentrate primarily on the "type II" model, 
where up-type quarks receive their mass from one Higgs doublet, and down- 
type quarks receive their mass from the second doublet. High energy constraints 
derived from the W mass, the full width of the Z and the bb partial width of 
the Z are combined with low energy constraints from T(b — > S7), T(b — ► ctv t ) 
and B°-B° mixing to determine the experimentally favored configurations of 
the model. This combination of observables rules out small charged Higgs 
masses and small values of tan/3, and provides some information about the 
neutral Higgs masses and the mixing angle a. In particular, constraints derived 
from the p parameter rule out configurations where the charged Higgs is much 
heavier or much lighter than the neutral Higgses. The agreement of the model 
with experiment is roughly as good as that found for the minimal standard 
model. We discuss a scenario where T(Z — > bb) is enhanced relative to the 
standard model result, which unfortunately is on the verge of being ruled out 
by the combination of T(b — > 57) and p parameter constraints. Implications for 
various extensions of the standard model are briefly discussed. 



I Introduction 



Precise measurements of electroweak observables are beginning to place significant 
constraints on extensions of the standard model. Measurements of the W mass 
the partial widths and forward-backward asymmetries of the Z |2[], and the branching 
ratio for B — * X s 7 [[3] all provide significant information concerning unknown pa- 
rameters of the standard model as well as constraining various types of new physics. 
Furthermore, recent experimental evidence for a top quark mass in the range 174 ± 17 
GeV [U eliminates much of the uncertainty in these constraints. To date, the minimal 
standard model has succeeded in its description of virtually all phenomena of elec- 
troweak origin || |J. At the same time, however, there are some sectors of the theory 
about which we know very little. In particular, predictions of the standard model 
depend very weakly on the mass of the Higgs boson, and as a consequence we can say 
very little about the Higgs sector of the theory. In part for this reason, the extension 
of the standard model to include a second Higgs doublet is at present far from being 
ruled out experimentally. Furthermore, the inclusion of a second Higgs doublet is a 
common feature of many extensions of the standard model, such as supersymmetry 
or certain versions of the SO(10) and SU(5) grand unified theories. It is therefore of 
interest to determine whether such a simple extension of the standard model Higgs 
sector is compatible with experiment, and what further extensions of the model are 
indicated if it is not. In addition, we would like to determine whether the discovery 
of a charged or pseudoscalar Higgs is necessarily an indication of supersymmetry, or 
whether such a discovery could simply be an indication of an extended Higgs sector. 

At tree level the two Higgs doublet model is identical, in most respects, to the 
standard model. It reproduces the important relation = Mf cos 2 ^, and fur- 
thermore the Yukawa couplings of the quarks can be chosen so as to eliminate flavor 
changing neutral Higgs interactions 0. The two Higgs model is also in agreement 
with experiment at the level of radiative corrections. The two Higgs model differs 
from the standard model, however, in that radiative corrections often depend rather 
sensitively on the details of the Higgs sector. For this reason, experimental data are 
beginning to significantly constrain the parameters of the Higgs sector of the two 
doublet model. Low energy data such as Y{b — > 57), T(b — > ctu t ) and the B° — B° 
mixing amplitude rule out small values of the charged Higgs mass, and significantly 
constrain tan f3, the ratio of the vacuum expectation values. High energy data such 
as the W mass and the full width of the Z constrain the Higgs masses and mixing 
angles; in particular, configurations which give a positive contribution to the p pa- 
rameter are disfavored by these data. Additional stringent constraints are provided 
by Rb = T(Z — > bb)/T(Z — > hadrons). Our aim here is to combine these constraints 
and to discuss the implications of a moderately heavy top quark for the model. 

In Sec. II, we present a general overview of the two Higgs doublet model, with 
particular emphasis on the structure of the Yukawa couplings. In Sec. Ill, we discuss 
oblique corrections in the model, and discuss their implications concerning the top 
quark mass. In Sec. IV, we introduce constraints on the model from high energy 
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data such as Rb and the full Z width. As a byproduct, we demonstrate that the 
combination of p parameter constraints with constraints from makes the model 
incompatible with a top mass larger than about 200 GeV. In Sec. V, we introduce 
constraints from low energy data such as V(b — > sj), r(b — > cru T ) and B° — B° mixing, 
and illustrate how these rule out small charged Higgs masses and small values of tan j3. 
In Sec. VI we combine all of the various constraints to determine the configurations 
of the model that are compatible with experiment for a top mass of 174 GeV. In 
Sec. VII we conclude. 



II The two Higgs doublet model 

In this section we briefly review some generalities of the two Higgs doublet model; for 
a more complete treatment, we refer the reader to Refs. 0, The two Higgs doublet 
model is essentially the SU(2) x U(l) standard model, modified by the extension of the 
Higgs sector to include two scalar isodoublets $i, $2, which we write in component 
form as 

* 4= ((# + fr?)/>/0" (1) 
The Yukawa couplings of the scalars and fermions are constrained by the requirement 
that there be no flavor changing neutral Higgs interactions. It was shown in Ref. 0] 
that a necessary and sufficient condition for the elimination of flavor changing neutral 
interactions in the Higgs sector is that each quark of a given charge must receive 
its mass from at most one Higgs field. One finds then that there are two possible 
arrangements for the Yukawa couplings: one may either couple all of the fermions 
to a single Higgs doublet (giving what is known as the "type I" model), or one may 
couple $1 to the right-handed down-type quarks, and $2 to the right-handed up- 
type quarks (giving what is known as the "type II" model). Here we will consider the 
latter possibility. Imposing the discrete symmetry $ 2 —* — $2, (u,c,t)n — > —(u,c,t) R 
is enough to ensure that the Yukawa couplings have the correct form. The Yukawa 
Lagrangian has the form 

33 3 

L Y = Y.kL^i,R + D i3%L^ldj,R + J2 U i3%L%Uj >R , (2) 

where £^1, qi t L are the left-handed lepton and quark doublets, U and D are the quark 
mass matrices, and $2 = —^2^*2- The Higgs potential, as in the standard model, has 
its minimum at non-zero values of the fields $1, $2- After applying the requirement 
that there be no flavor changing neutral interactions, one finds that the vacuum 
expectation values (VEVs) can be chosen both real and positive: 



( $ i)= .,.,./o > <*»> = .,_,"/o • (3) 
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Table I: Quark Yukawa couplings in the "type II" two doublet model. Overall factors 
of —ig/2M w , —g/2M w , igV ud jM w have been omitted in the couplings of Hf 2 , if 3, 
and H + respectively. 



Higgs 


uu 


dd 


ud 


Hi 
HI 

HI 
H+ 


m u sin a / sin (3 
m u cos a/ sin (3 
m u cot (3% 


rrid cos a / cos (3 
—m d sin a/ cos (3 
m d tan/?7 5 


m u cot/3(^)+m d tan/3(^f) 



After diagonalizing the Higgs mass matrix, one finds that the fields <pf mix to form a 
charged Nambu-Goldstone boson G + and a charged physical scalar H + of mass m + : 

G + = cos^ + sin/30+ F + = -sin/30+ + cos/3^. (4) 

Similarly, the imaginary parts of the neutral components Xi 2 m i x to form a neutral 
Nambu-Goldstone boson G° and a CP-odd physical scalar Hi, again with mixing 
angle j3: 

G° = cos/5 X ? + sin/? X °, P 3 = -sin/5 X ? + cos/?X2- (5) 

Finally, the scalars (f)\ 2 m ix to form a pair of neutral CP-even scalars, now with 
mixing angle a: 

H^ = cos a4>\ + sin acj)^ H% = — sin + cos ac^. (6) 

The masses m 12 of obey m\ > m 2 . 
The mixing angle /3 is given simply by 

tan/5 = — . (7) 

Vl 

The quantity tan (3 plays a central role in the theory because the Yukawa couplings 
are often proportional to either tan (3 or cot j3. The Yukawa couplings of the various 
quarks to the various scalars have been summarized in Table |. We see in partic- 
ular that couplings of up-type quarks are enhanced for small values of tan /3, while 
couplings of down-type quarks are enhanced for large values of tan (3. Furthermore, 
one can decouple Hi from the up-type quarks by choosing a = 0, and so on. The 
large top quark mass, in combination with the enhancement of the top Yukawa cou- 
pling for small values of tan/3 makes it possible to derive a lower bound on tan/3. 
Furthermore, we note that since up-type quark masses are proportional to f 2 and 
down-type masses to V\, the mass hierarchy mt 3> rrib tends to favor large tan/3. 
In the type I model, where all of the quarks receive their masses from (say) $1, the 
Yukawa couplings of the fermions are all proportional to cot (3 or esc /3; consequently 
the Yukawa couplings of the bottom quark are negligible compared to those of top 
quark regardless of the value of tan j3. 
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Ill Oblique Corrections in the Two Higgs Dou- 
blet Model 



In this section we discuss some features of oblique corrections to electroweak observ- 
ables in the two Higgs doublet model. 

It has been known for quite some time that oblique radiative corrections from the 
Higgs sector can be significantly larger in the two Higgs model than in the minimal 
standard model || [10], In particular, the contribution of the Higgs sector to Ap, 
defined by 

EwwiO) £zz(0) , , 

a p = — ' w 



Ml 



where ^ww,zz are the W and Z self-energies, can be significantly larger in the two 
doublet model than in the standard model. The non-standard contributions to the 
vector boson self energies in the two doublet model have been presented in Refs. H, 



12], and are summarized in the Appendix. In the minimal standard model, the 



contribution of the Higgs to Ap is given, in the limit mm 



oo, by 



ApHiggs (MSM) 



3G P M 2 



w 



tan 2 6 W 



log- 



m 



Higgs 



(9) 



and grows only logarithmically with the Higgs mass, 
contribution to Ap is given by 



By contrast, the top quark 



Ap t0 p 



3G, 



-m 



t ■ 



(10) 



which grows quadratically with m t . It is the weak dependence of Ap on m^ggs that 
excludes top quark masses larger than 200 GeV in the minimal standard model. Since 
Ap grows quadratically with m t , but falls off only logarithmically with mHig gs , one 
must have an exponentially large Higgs mass to prevent Ap from becoming too large 
when the top is heavy; since we expect on general grounds that m H igg S < 1 TeV, the 
standard model cannot accommodate a top mass larger than about 200 GeV. In the 
two Higgs doublet model, the situation is different. Here the leading behavior (for 
Higgs masses much larger than the W mass) of Ap is given by ||TT] 



Ap 



3G, 



Higgs 



8v^7T 2 



sin 2 (a — P)F(m + , m 3 , mi) + cos 2 (a — (3)F{m +1 m 3 , m 2 )^ , (11) 



where 



F(m + ,TO3,TOi i2 ) 



m , — 



2 2 

m + — mj 



m.m 



+ '"1,2 



m\ — m\ 2 



log 



i ( m + 

2 



?77 + 

m? 2 



+ 



mj - m\ 2 



log 



m 3 

mj 2 , 



(12) 
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The Higgs contribution to Ap can be either positive or negative depending on the 
ordering of the Higgs masses, and in general grows quadratically with the largest 
Higgs mass. If the Higgs masses are ordered as 

< m + < m 3 or m 3 < m + < m^, (13) 

then ApHiggs will be negative and grow quadratically with the largest Higgs mass; if 
the charged Higgs is heavier or lighter than all of the neutral Higgses, then ApHi ggs 
will be positive. The negative contribution to Ap is largest if m + ~ 0.562mh ea vy, 
where m heavy is the largest Higgs mass. For appropriate values of a — ft, one then has 

Apmggs ^ x 0.216mL vr (14) 

Since the top quark contribution to Ap is quite large, the Higgs contribution to Ap 
must be either small (if it is positive) or negative. 



IV Constraints from the full Z width and 

In this section we introduce constraints derived from the full Z width and the bb 
partial width. The LEP measurements of these quantities are 

T(Z -> all) = 2.4974 ± 0.0038 GeV (15) 

and 

T(Z 66) 

^(z Uadrlr 2202 ±0 002 °- <16) 
The quantity Rb is particularly convenient to work with, since oblique and QCD 
corrections to T(Z — > bb) and YiZ — > hadrons) cancel to a large extent in the ratio R^. 
On the other hand, the full Z width is rather sensitive to oblique corrections through 
the p parameter. We begin with a general discussion of the radiative corrections to 

nz - /./•). 

The Z width in the two Higgs doublet model has been studied previously in 
Ref . , and, in the context of the supersymmetric standard model, in Ref . (14| . 



We have independently carried out the calculation, and our results agree with those 
found previously. The vertex connecting the Z to a pair of fermions may be written 

as 

T tt = i&y ll {v -075), (17) 
where at tree level the vector and axial couplings v and a are 

v = and a = . (18) 

2sc 2sc v ; 

We use the abbreviations s = sm9 w , c = cos^vf- Radiative corrections modify these 
couplings through self energy and vertex corrections. The self energy corrections arise 
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from the Z self energy and from the Z-photon mixing. The Z self energy insertion 
introduces an overall factor of the Z wavefunction renormalization constant. The 
Z-photon mixing modifies the vector coupling of the Z to a fermion pair: we make 
the replacement 

h-2Qs 2 , ^ AZ (M 2 ) 



2sc 



+ Q- 



(19) 



to include this effect. Vertex diagrams introduce further corrections to the vector and 
axial couplings of the Z, which we denote by 5v, 5a. Neglecting fermion masses, we 
find that the partial Z width into a particular final state may then be written 



T(Z - //) = Z z 



n c aM z 



v 2 + a 2 + 2v Re 5v + 2a Re 5a + 2vQ Re 



Xaz(M 2 ] 
M 2 



(20) 



where n c is the number of colors and Zz is the Z wavefunction renormalization 
constant, which cancels in the ratio Rb- We have employed the renormalization scheme 



of Hollik |y| throughout our calculations. 

In addition to these purely electroweak corrections, partial widths into quarks are 
modified by QCD corrections. For quarks other than the b, the QCD corrections 
enhance the partial width by about 4%: we have 



T(Z - qq) = T°(Z - qq) 



1 + 



a s (M : 



+ 1.41 



7T 



a s (M : 



7T 



(21) 



where T is the uncorrected width. In the case of b quarks, the finite b mass and the 
large mass of the top quark modify the QCD corrections. The relevant corrections 
have been calculated in Refs. [T7|, [RJ. In Ref. it was pointed out that the 
vector and axial parts of the Z width into b quarks are affected differently by QCD. In 
particular, the vector-current contribution to the Z width into b quarks is corrected 
by the same QCD factor appearing in Eq. (|2T|), while the axial part Ta is corrected 
by a factor 



1 A 



1 + 



7T 



H 1.41 - J( 



7T 



(22) 



where z = M^/Am 2 , and —I(z)/3 is negative and increases in magnitude for increas- 
ing top mass. For a top mass of 100 GeV, —I(z)/3 ~ —3; for 250 GeV, —I(z)/3 ~ —5. 
The impact of these QCD corrections is to slightly reduce Rb- A second set of QCD 

T7| , [T5[| , where corrections to T(Z — > bb) due 
These corrections enhance YiZ — > bb) by 
. The combined effect of all of these QCD 



corrections has been discussed in Refs 
to the finite b mass have been calculated, 
about 2 MeV, giving a small increase in Rb 



corrections is to enhance Rb by about 0.4%. 

In Fig. 1 we have plotted Rb as a function of the top mass for both the minimal 
standard model and the two Higgs model for a few values of tan j3. 

The Z width into b quarks is of particular interest because of the presence of 
virtual top quarks in various vertex corrections. These corrections can be used to 
constrain tan/? and the charged Higgs mass, and, to a lesser extent, the mixing angle 
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Figure 1: ^ as a function of the top mass in the minimal standard model 
and in the two Higgs model. The solid curve is the standard model result for 
m ffiggs — 100 GeV. The upper and lower dashed curves are the two Higgs results 
for m 2 = m 3 = 50 GeV, mi = 875 GeV, m + = 422 GeV and tan /3 — 70, 1 respec- 
tively. For comparison, the standard model result for R d , the corresponding quantity 
for d quarks, is also shown. The error bars indicate the la experimental measurements 
of m t and Rj,. 
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Figure 2: Vertex diagrams involving charged and neutral Higgs exchange in the two 
doublet model. Diagrams involving neutral Higgs exchange that do not contain po- 
tentially large trigonometric factors (such as cot 2 f3, sec 2 0) have been neglected. 
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a and the neutral Higgs masses. The Higgs exchange vertex diagrams contributing to 
YiZ — > bb) are shown in Fig. 2. These break down into two classes: those involving 
charged Higgs exchange, and those involving neutral Higgs exchange. The dependence 
of Rb on the top mass results from the diagrams involving charged Higgs exchange, 
since these involve the top quark Yukawa couplings of Table |I[ These diagrams tend to 
suppress the decay of the Z into b quark pairs by an amount that grows quadratically 
with the top mass. In the minimal standard model, the effect of a heavy top is to 
suppress the Z width into b quarks by an amount [19[ 



8ns 2 As 2 c 2 V 3 



+ U + log 



V3 6cV 



(23) 



in the limit m t 3> My/- In the two Higgs model, diagrams involving charged Higgs 
exchange further suppress the Z width into b quarks by an amount (for m + , m t ^> 
M w ) 

Arty ^ a QEpMz 1 m 2 ( 2s 2 \f x x 
AY(Z^bb) H ± ~ o „ — , — 2 „ a 1-— --— — i logx + 



g7rs 2 4s 2 c 2 M 2^ tan 2 £ y 3 (X - l) 2 ° X - 1 

where x = m 2 /m 2 + . This correction falls off for large m + and grows with increasing 
m t . The corrections due to charged Higgs exchange are always negative, and increase 
with increasing top mass. This makes it difficult to accommodate a 150-200 GeV top 
quark mass, particularly for small values of tan /3; for a top quark mass of 174 GeV, 
these diagrams exclude small values of tan f3. 

The situation is somewhat different in the case of diagrams involving neutral Higgs 
exchange. Here we find that the contribution to YiZ bb) can be either positive 
or negative, depending on the masses of the neutral scalars. These diagrams have 
no top mass dependence, of course, since they involve only 6-quarks. Their explicit 



forms have been reported in Ref. JT3[, and are summarized in the Appendix. These 
diagrams are only important for large values of tan/3, for it is in this situation that 
the bottom quark's Yukawa coupling becomes large. In Fig. 3 we have plotted the 
correction to YiZ — > bb) resulting from neutral Higgs exchange as a function of 
for the special case a = tt/2, rri2 = 50 GeV, and tan/3 = 70. In this case H® 
decouples from the vertex. The correction is positive for small and roughly equal 
masses m<i ~ ms, but becomes negative for large mass splittings. This is in fact 
one of the few situations where vertex corrections significantly enhance the Z partial 
width into b quarks. These positive corrections make it possible to accommodate a 
heavy top quark, particularly if the neutral Higgs bosons are light and tan (3 is large. 
As we will show below, however, large top quark masses are not admissible in the 
model. 

We have also studied the impact of these vertex corrections on the partial widths 
for Z — > t + t~ and Z — > v T v T . Here the Higgs exchange vertex corrections are 
important only for large values of tan/3 and when two or more neutral scalars are 
light. Both of the partial widths are enhanced by a small amount by these corrections; 
however, the effect is less than 1 MeV in both cases. 
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In light of the present discrepancy between the standard model prediction of Rb 
and the measured value, one interesting feature of the two doublet model is the 
enhancement of Rj, relative to the standard model result when tan/? is large and 
two or more of the neutral scalars are light. Although we will show below that such 
configurations of the model are on the verge of being ruled out by F(b — > sj), we 
feel that direct experimental tests of this scenario are worth briefly discussing. Since 
the enhancement of Rj, requires small masses for the neutral scalars, it may be rather 
easy to test at LEP II. The cross section for the process e + e~ — > Z*H^ is suppressed 
relative to the standard model result for e + e~ — > Z*H° by a factor of sin 2 (a — (3). If 
sin 2 (a — (3) is not too small, it should be possible to detect a light scalar in the mass 
range rri2 ~ 50 — 100 GeV at LEP II using data from a standard model Higgs search. 
Barring extremely small values of sin 2 (a — /3), this should provide a fairly conclusive 
experimental test of this scenario. In the event that sin 2 (a — (3) does turn out to be 
small, the process e + e~ — > Z* — > H^H®, for which the cross section is proportional 
to cos 2 (a — (3), may provide a complementary experimental probe of this scenario. 
The current lower bound on m<i |20| ranges from about 58.4 GeV for sin 2 (a — 0) — 1 
to about 30 GeV for sin 2 (a — (3) = 0.05. Finally, we note that this scenario is 
relevant only to the type II model. In a type I model, the Yukawa couplings of both 
top and bottom quarks are proportional to cot (3. Consequently the positive vertex 
corrections from neutral Higgs exchange are far smaller than the negative corrections 
from charged Higgs exchange, and so one does not expect significant enhancements 
of YiZ — > bb) at large values of tan/3. 

The various corrections to the Zbb vertex also modify the bb forward-backward 
asymmetry of the Z. At tree level, Afb(bb) is given by 

A fUU\ 3 2V ^ e 2VbCLb fOK\ 

A ' m = I vl + alvl + aV (25) 

where Vf, a/ are the vector and axial couplings of the electron and b quark. Vertex 
corrections tend to modify the left-handed coupling of the b for small values of tan/3, 
and the right-handed coupling of the b for large values of tan (3. We have calculated 
the bb forward-backward asymmetry of the Z including oblique and vertex corrections; 
additional contributions from box diagrams are known to be numerically small and so 
are neglected here. The forward-backward asymmetry shows much the same behavior 
as Rb as a function of the Higgs masses and mixing angles: Af b (bb) is suppressed at 
small values of tan/3 as a result of charged Higgs effects, while in the case of large 
tan j3 there are several possibilities. If all of the neutral Higgses are light (with masses 
less than 100 GeV) and tan/3 is large, then Afb is enhanced slightly; if H®,!!® are 
light, but H® is much heavier, then Afb is enhanced for \a\ ~ ir/2, but suppressed 
for a ~ 0. In addition to its sensitivity to vertex corrections, Afb is also sensitively 
dependent on oblique corrections through the value of sin 2 9w The enhancement of 
Afb at large values of tan/3 limits the extent to which R b can be enhanced by neutral 
Higgs effects. 

To illustrate the range of top masses that can be accommodated by the two doublet 
model, we can attempt to construct a scenario which allows a very large top mass by 
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Figure 3: Correction to T(Z — > bb) from exchange as a function of m 3 for 
m 2 = 50 GeV, tan/3 = 70, and a = tt/2. 
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arranging for negative contributions to Ap and positive contributions to Rb from the 
Higgs sector. In order to have positive contributions to Rb, m 3 must be small. The 
requirement of a small mass for together with the fact that m 2 < mi implies that 
we must have the mass hierarchy m 2 ,m 3 < m + < m\. We will then have a negative 
contribution to Ap for \a — (3\ = it/2. Here, however, we run into a conflict with the 
the constraint from R^. in order to have simultaneously tan/3 large and \a — (3\ = n/2, 
we must have (3 ~ n/2, and a ~ 0. Referring to Table [I], we see that for this value of a 
the Yukawa coupling of is small; hence it is not possible to simultaneously obtain 
a negative contribution to the p parameter and still have positive vertex corrections 
to Rb- Although either constraint separately can permit a very large top mass (on the 
order of 250 GeV or larger), when the two are combined the top mass is constrained 
to be less than about 200 GeV, consistent with both standard model expectations 
and recent evidence for top quark production at CDF H. The existence of this upper 



bound has been noted previously by the authors of Ref. [^TJ. The mass limit of 200 
GeV has been established by a systematic search of the entire parameter space of the 
two doublet model. 

The full Z width provides constraints on the charged Higgs mass. The Z width 
is quite sensitive to oblique corrections through the p parameter: if the Z width is 
expressed in terms of Gf, the effect of oblique corrections is to shift the value of sin 2 0w 
and to renormalize the width by a factor p. As a result, positive contributions to Ap 
tend to enhance the Z width. This constrains the mass splittings in Higgs sector of 
the two doublet model; in particular, at the la level and with a 174 GeV top quark, 
we find that the charged Higgs mass is constrained by 

flight — 130 GeV < m + < m heavy + 130 GeV, (26) 

where m heavy light are the masses of the heaviest and lightest neutral Higgses, respec- 
tively. These mass limits have again been derived by a systematic search of the entire 
parameter space of the two doublet model, allowing each of the Higgs masses to vary 
between 50 and 1000 GeV, and each of the mixing angles to vary over their full range. 



V Constraints from Y(b —> S7), B —B mixing, and 

r(6 — > ctv t ) 

The CLEO collaboration || has recently reported a value for the inclusive branching 
ratio for radiative B decays: BR(5 -> X s -f) = (2.32±0.51±0.29±0.32) x 10^ 4 , where 
the first error is statistical, and the latter two arise from systematic errors in the yield 
and efficiency, respectively. The 95% confidence level limits on the branching ratio 
are 

1 x 10" 4 < BR(B -> X s7 ) < 4 x 10" 4 . (27) 

This result can be used to rule out small charged Higgs masses and small values 
of tan/3. The relevant electroweak diagrams are shown in Fig. 4, and have been 
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calculated in Ref . |22|] . Ignoring QCD corrections, the effective Lagrangian giving rise 
to the b — > 57 transition is 



G F 

Lett = _ ~^j2^ th ^~ 



ts o 2 



Ml, 



2 



(28) 



where a is a color index, and F, G are given by 

3x 3 



F(x) 

G{y) 



4(x — I) 4 
3y 2 - 2y 



2x 2 , -8x 3 - 5x 2 + 7a; 

logx + 



6(y - 
+ cot 2 /? 



1)3 

3y 3 



logy + 



24(x - l) 3 
-5y 2 + 3y 



2y 2 



12(1/ -1)' 



12(y - ] 
logy + 



-8y 3 - 5y 2 + 7y 



72(y - If 



(29) 



These results are only valid when one ignores QCD effects. In order to incorpo- 
rate QCD, it is necessary to resum the large logarithms ~ log(M^/m 2 ) using the 
renormalization group. The leading logarithmic QCD corrections result in an overall 
enhancement of the b — > 57 transition rate by a factor of 3 to 5 over that which 
would be obtained using this naive effective Lagrangian. These corrections have been 



discussed extensively in the literature [|23| , [24] , P5| , |26| , p7| . Here we will use the QCD 
calculations of Ref. |[25|| . There are minor differences between the various calculations 
of the anomalous dimension matrix; however these result in only small (on the order 
of 1% or less) variations in the calculated value of BR(6 — » 57). To extract a value of 
BR(i9 — > X s j), it is convenient to use the approximation f23| 



T(B -> X s7 ) r(6 -> 57) 



T(B — > X c eh> e ) Y(b — » cei/ e 



(30) 



which reduces uncertainties due to a factor m| which appears in both numerator and 
denominator. Including the leading logarithmic QCD corrections, we then have 



T(B -> X s7 ) 
r(S -► Xez/ e ) 



6a 


v; s v tb 


2 


iif{m c /m b ) 




2 



\A{m b )\< 



(31) 



where A(m&) is the coefficient of the operator (e/87T 2 )m b s lct ^b^F^ evaluated at the 
scale m b and f(m c /m b ) ~ 0.316 is the phase space factor for the semi-leptonic decay. 
The coefficient A(m b ) is given in terms of r\ = a s {m b ) / a s (M w ) by 



A(m b ) = ?7"23A(Mi 



w. 



8 . 



14 

23 



,1 j 



(32) 



where A(Mw) = F(x) + G(y), B(Mw) is the coefficient of the "chromo-magnetic 
moment operator" (g/87T 2 )m b s1a, lu T^ l3 b l lG^ u evaluated at the scale Mw, and the 
sum results from mixing with the four quark operator {cLal^bi^isLpl^CLp)- The 
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Figure 4: Diagrams in the two Higgs model contributing to the magnetic moment 
type operator that mediates the b — *> 37 transition. 





numbers aij and Xi are related in a simple way to the eigenvalues and eigenvectors of 
the anomalous dimension matrix. We have used | V^V^,| 2 /| Vb c \ 2 = 0.95. 

Recent analyses of the uncertainties in the leading logarithmic calculation have 
been given in Refs. |27], 28[1 , where it was pointed out that the theoretical uncertainty is 
dominated by the unknown next-to-leading logarithmic corrections. Other significant 
uncertainties result from the experimental error in the measurements of a s [Mz) = 
0.12 ±0.01 and the ratio m c /mfc, which occurs in the phase space factor for the semi- 



leptonic decay. Combining the various theoretical uncertainties given in Ref. [27 



in quadrature, one finds that the overall error for the standard model result is on 



the order of 30%; the error estimate of Ref. is also roughly of this magnitude. 
In the case of the two doublet model, the situation is complicated by the greater 
sensitivity of the result to the next-to-leading order corrections. This has been 
pointed out in Ref. [ 2"8fl . In light of this, the constraints from b — > 57 presented here 



should be considered as qualitatively correct, but nonetheless subject to considerable 
theoretical uncertainty. To factor in these uncertainties, in deriving limits on tan/3 
and the charged Higgs, we have employed a 30% estimate of the theoretical error. We 
require only that 

0.7 x (Theoretical Estimate) < Experimental Upper Bound. (33) 

This has the effect of weakening somewhat the lower bound on the charged Higgs 
mass. 

Additional low-energy constraints on the two doublet model can be obtained using 
B° — B° mixing. Here the relevant quantity is Xb = Am/T, where Am is the mass 
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difference between the heavy and light admixtures of B° and B°, and T is their average 
width. By observing time-dependent B° — B° oscillations, ALEPH and DELPHI 
have determined the mass splitting to be Am = 3.41 x 10~ 4 eV, which yields Xb = 
0.797±0.128 [ 29[ . Measurements of the time integrated mixing probability by ARGUS 
and CLEO give Xb = 0.67 ± 0.08 p9f. The weighted average of these values is 



Xb = 0.706 ± 0.068. In the two doublet model, Xb is given by 



Am G\ 



x h 



where r B o 



T 67r 2 
1.53 ±0.09 ps 



/? %* d \ 2 \Vt b \ 2 f B m B B B rj B T B M^ { I ww + I WH +I HH ), (34) 



is the B lifetime, f B is the B meson decay constant, 
m B is the B meson mass, B B is the bag factor, and r\ B is a QCD correction factor, 

The CKM matrix element 



whose value lies in the range from 0.55 to 0.85 5T[ |32 
\V t d\ lies in the range from 0.005 to 0.014 |B3], while V t 



tb 



on m+, m t and tan/3 is contained in Iww, Iwh, and I B h 



1. Finally, the dependence 
, We have HI 



iww 



' WH 



' HH 



1 + 



3 — 9x 6x 2 logx 



xy cot 2 P 
xy cot 4 P 



X- l) 2 

(4z 



+ 



(x-1) 
1) logy 



3 logx 



4 



-z) 2(1 
2y log |/ ' 



where x 



2(l-y)2(l 

ll ' d //! : 

m 2 /mi, and 2 



+ 



x - 4 



x) 



l-z) (l-x)(l-y) 



and 



(35) 



Ml/ml 



m t/Mw, y — "hi aiLu - * — ±y± wi 

A final constraint on the charged Higgs mass and tan P is provided by T(b — > ctv t ). 
Since in the two doublet model this process is mediated by both W and charged Higgs 
exchange, one cannot have too small a mass for the charged Higgs or too large a value 
of tan P without exceeding the experimental limits. This constraint has been analyzed 
recently in Ref. [Q, where it was shown that the allowed values of m + and tan P must 
satisfy 

tan/3 < 0.51m+ [GeV] (36) 

at the la level. In contrast to the constraints given above, this measurement provides 
an upper bound on tan P, rather than a lower bound. 

In Fig. 5 we have plotted the allowed region in the m + -tan/? plane after the 
constraints due to BR(6 — > sj), T(b — > cru T ), B° — B° mixing, and Rb have been 
applied, with a top mass of 174 GeV. The constraints from Rb indicated are at the 
2.5cr level; to agree with experiment to within 2a or less, tan/? must be large and 
the scalars H® 3 must be light. In this situation, the charged Higgs mass cannot be 
much larger than about 150 GeV, for reasons we will discuss in Sec. VI. We have used 
the 95% confidence limits for b — > s'j. The constraint from BR(6 — ► 57) rules out 
small charged Higgs masses, while the inclusion of Rb strengthens the lower bound 
on tan p. The constraint from Xb is similar to that from Rb but somewhat weaker. 
Similar conclusions have been reached in Refs. [571 EBL From these considerations, 
we can constrain m + and tan P by 



m+ > 200 GeV 



(37) 
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and 



tan/3 > 0.7 



(3* 



at the 2.5cr level. At the 2a level, the bound on tan j3 is roughly 13 for m t = 157 GeV, 
increasing to 38 for m t = 174 GeV. 

Since experimental measurements of BR(6 — > 57) are likely to improve, we present 
also the charged Higgs mass limits that result for smaller values of BR(6 — > sj). This 
serves to illustrate the high sensitivity of these mass bounds to theoretical uncertain- 
ties. For comparison, we give both the mass bound derived using the leading order 
QCD corrections we have employed here, as well as the mass bounds reported in 
Ref. [28[ where the next-to-leading order corrections were partially included: 



f 200 GeV (LO), 150 GeV (NLO) BR(6 -> 57) < 4 x 10" 



m + > { 380 GeV (LO), 200 GeV (NLO) BR(6 -> 57) < 3 x 10' 



l 650 GeV (LO), 300 GeV (NLO) BR(6 -»s 7 )<2x 10" 



(39) 



We see that the lower bound on the charged Higgs mass may move downward by a 
factor of two once the full next-to-leading QCD corrections are known. We see also 
that the partial inclusion of next-to-leading order corrections relaxes the charged 
Higgs mass bound from the 200 GeV cited above to about 150 GeV. 



VI Global constraints on the two doublet model 

In this section we discuss the experimentally preferred configurations of the two dou- 
blet model by combining the constraints of the previous sections. Throughout we 
assume a top quark mass of 174 GeV. We first discuss those configurations of the 
model which are compatible with each of the various high energy measurements 
T(Z — > all), and M w ) at the 1.5a level. We then add to these the low energy con- 
straints. 

To be consistent with the LEP measurement of Rb at the 1.5<r level, one must 
assume that tan/3 is large and that and H® are quite light, having masses less 
than or of the order of 100 GeV. We have evaluated the maximum allowed masses 
and the minimum value of tan (3 by both a systematic search of the parameter space 
and by randomly sampling a large number of configurations of the model. By both 
methods, we have found the limits 

m 2 ,m 3 < 65 GeV, (40) 

and 

tan/3>53, (41) 

which translates to a value of (3 very nearly equal to tt/2. In addition, provides a 
constraint on the mixing angle a in the case that m\ is large; the reason is that, for 
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Figure 5: Allowed region in the m + -tan/3 plane for the combination of observables 
indicated. The shaded region is allowed. 
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mi large and a ~ 0, the neutral Higgs contribution to Rb is negative (cf. Fig. 3). As 
a consequence, one has 

\a\ > vr/3 for m x > 100 GeV. (42) 

If one assumes the values of m 2 3 and tan f3 implied by Rb, then further constraints 
can be derived from M\y and the full Z width, through the p parameter. The leading 
behavior of Ap in this case is given by 

Apm - * {< ~ cos2 log G|) ) ' (43) 

From the leading behavior of ApHi ggs , we see immediately that the charged Higgs 
mass cannot be too large, for this would give a large positive contribution to Ap. 
Furthermore, although the second term in Eq. is negative, its coefficient is small 
due to the bound on a imposed by Rb when mi is large. As a result, we have 

m+ < 150 GeV. (44) 

Due to the quadratic growth of Ap with m + , larger values are strongly disfavored. 
Finally, large values of mi are weakly favored by p parameter constraints; one has 
mi > 500 GeV, although smaller values are not in terribly poor agreement with the 
data. In combination with constraints from Rb, the bound on the mass of mi gives 

|a| > tt/3. (45) 

This configuration of the model is also compatible with the bb forward-backward 
asymmetry of the Z. In deriving these limits, we have again carried out a systematic 
search of the parameter space of the model, covering the full range of a and j3, and 
varying each of the scalar masses between 50 and 1000 GeV. If one allows larger 
values of m 1; then slightly larger values of m + are admissible; for m x = 1.5 TeV, the 
upper limit on the charged Higgs mass recedes to 170 GeV. 

Finally, we can combine the constraints from high energy data with those from 
b — *■ sj, b — > crv T , and B mixing. This class of constraints also rules out small values 
of tan/3, although the bounds are weaker than those derived from Rb. If we take 
the weaker of the two charged Higgs mass bounds given in Eq. (^), we find that the 
charged Higgs mass is constrained to lie in a narrow region in the vicinity of 150 GeV: 

m + ~ 150 GeV. (46) 

We find that compatibility of the model with the measured value of Rb, when com- 
bined with the other relevant observables, yields a very restrictive set of constraints. 



VII CONCLUSIONS 

We have presented constraints on the type II two Higgs doublet model from a variety 
of observables. The overall agreement of the model with experiment is acceptable, 
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and quite comparable to that found for the minimal standard model. We have shown 
that tan /3 and the charged Higgs mass are bounded by 

m + > 150 ~ 200 GeV, tan/3 > 0.7 (47) 

for 2.5cr agreement with measurements of R^, and discussed the sensitivity of the 
charged Higgs mass bound to theoretical uncertainties. The lower bound on tan/3 
increases to about 13 for 2cr agreement with Rj,. Oblique constraints on the model 
require that the charged Higgs mass lie in the range 

mught - 130 GeV < m+ < m heavy + 130 GeV, (48) 

where m heavy light are the masses of the heaviest and lightest neutral Higgses. We have 
discussed one very tightly constrained scenario where the model (unlike the standard 
model) agrees with measurements of R^ at the 1 to 1.5<r level. In this case, every 
parameter of the model is constrained in some way; for the mixing angles one has 

tan/3>53, \a\ > tt/3, (49) 

while for the Higgs masses one finds 

m 2j3 < 65 GeV, m + ~ 150 GeV, m 1 > 500 GeV. (50) 

The preference for large values of m x is weak, and p parameter constraints strongly 
disfavor larger values of m + . This configuration of the model is virtually ruled out 
by b — > sj: only a very narrow allowed region remains for the charged Higgs mass, 
and if the uncertainty on the measurement is reduced even slightly, this window is 
likely to disappear. Consequently, it is improbable that the two doublet model can 
improve on the standard model prediction of R^. 

The large tan/3 scenario discussed here may nonetheless be of some interest for 
studies of supersymmetric grand unified theories. Indeed, if one assumes a large 
value of tan/3 and a charged Higgs mass on the order of 100 GeV, then one finds 
that i^23 are light and nearly degenerate, and that a is very nearly equal to —tt/2. 
This configuration of the model is surprisingly similar to that found in the large tan /3 
scenario. In supersymmetric SO (10) models the Yukawa couplings of the t, b, and r 
can be made to unify at the unification scale if tan/3 ~ 50 — 60 and m t ~ 160 — 190 
GeV |36|, ^]. Of course, the full calculation including superpartners is necessary in 
order to study the viability of this scenario. 

We have shown that the type II two Higgs doublet model is roughly consistent 
with numerous experimental observations at the level of radiative corrections. The 
measured values of the W mass, the Z width, and low energy data such as T(b — *> 57), 
r(6 — > ctu t ) and the B° — B° mixing amplitude can all be accommodated within the 
model. The partial width of the Z into b quarks poses a significant challenge to the 
model, which will become increasingly restrictive as data on this and other observables 
improve. We conclude, then, that the minimal two doublet model without supersym- 
metry cannot be ruled out on the basis of current data. Although the discovery of 
a light neutral Higgs or a charged Higgs would be indicative of supersymmetry, it is 
not an ironclad guarantee that the world is indeed supersymmetric. 
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IX APPENDIX 



In this appendix we collect the vector boson self energies and vertex corrections 
necessary for our calculations. The vector boson self energies can be expressed in 
terms of the two point integral 



167T 2 



I (p 2 , mi ,m 2 ) = —^ d D k 



1 



{k 2 — m 2 )([k — p] 2 — m 2 ) ' 



Another frequently appearing combination is 



AI (p 2 ,m 1 ,m 2 ) 



I (p 2 , m 1 , m 2 ) - J (0, mi, m 2 ) 

p2 



(51) 



(52) 



For the vertex corrections we need in addition the vector two point integral I±, defined 
by 

167T /<n\ f , n, k,. 



p ll I 1 (p l ,m 1 ,m 2 ) = ^-^ jd D k 



the scalar three point integral K , 



Kn = 



16tt 2 



d 4 k- 



{k 2 — m 2 )([k — p] 2 — m 2 ) ' 
1 



(k 2 -m 2 )([k-p 2 } 2 -m 2 )([k + Pl ] 2 -ml)' 
and the tensor three point integrals ^11,22,12,00 defined by 

PiijPivKu + p 2 ^p 2u K 22 + (p llM p 2u + PiuP2^)K 12 + g^Koo = 



(53) 



(54) 



f d D k 



kp}z v 



[k 2 — m 2 )([k — p 2 ] 2 — m 2 )([k + pi} 2 — m\ 



(55) 



The arguments of the three point integrals are (q 2 , p\, p\, mi, m 2 , m%) , where pi,2 are 
the momenta of the final state fermions and q = pi + p 2 . The arguments of these 
integrals are given below in the form (q 2 ,mi,m 2 ,m 3 ). The unrenormalized vector 
boson self energies S^a, ^az, ^zz, and Y* ww are 

T, AA (p ) = 



a 
An 



Z A z(p 2 ) 



{aim u r 
E — [(P 2 + 2mJ)/ (p 2 , m f , m f ) - 2mJ/ (0, m f , m f ) - p 2 /3 

(3p 2 + AM^)I (p 2 , M w , M w ) - 4M^/ (0, M w , M w ) 

(Am 2 + - p 2 )I (p 2 , m + , m+)/3 - Am 2 + I (0, m + , m+)/3 - 2p 2 /9 }, (56) 

AN f c v f Q f 



a 
Ait 
1 

3sc 



E 



(p + 2m f )I (p ,mf,mf) - 2m / J (0, m /; mj) - p /3 



[9c 2 + l/2]p 2 + [12c 2 + 4]M^ J (j9 2 , M w , M w ) 



12c 2 - 2)M^/ (0, M w , M w ) + p 2 /3 

(2m 2 + - p 2 /2)I (p 2 , m+, m+) - 2m 2 + I (0, m+, m+) - p 2 /3 



c 2 -s 2 



3sc 



>, (57) 
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Szz(p 2 ) 



4ttIY 3 



; / + a /) + 2m 2 f ]I (p 2 ,m f ,mf) ~ 2m 2 / (0, m f ,m f ) - p 2 /3 



m 



f 



2s 2 c 2 
1 

6s 2 c 2 



Io(p 2 ,m f ,m f ) 



[24c 4 + 16c 2 - 10]M^ + [18c 4 + 2c 2 - l/2]p 2 ) I (p 2 , M w , M w ) 



(24c 4 - 8c 2 + 2)M^/ (0, M w , M w ) + (4c 2 - l)p 2 /3 
1 

12s 2 c 2 



(c 2 — s 2 ) 2 ([4m + — p 2 ]Io(p 2 , m + , m + ) — Am+Io^O, m + , m + ) — 2p 2 /3^j 

+ cos 2 (a - (3) ([2ml ~ lOMf - p 2 ]/ (p 2 , m u M z ) - [ml - M 2 z ] 2 AI (p 2 , m u M z ) 

- 2m 2 / (0, mi, mi) - 2A/f J o (0, M z , M z ) - 2p 2 /3) 

+ sin 2 (a - 0) ([2m 2 - lOMf - p 2 }I (p 2 , m 2 , M z ) - [m 2 - Mf ] 2 A/ (p 2 , m 2 , M z ) 

- 2m 2 / (0, m 2 , ma) - 2M 2 / (0, M Z) M z ) - 2p 2 /3^ 

+ sin 2 (a; — /3) ([2m 2 + 2m 2 — p 2 )Io(p 2 , mi, m 3 ) — [m 2 — m 2 ] 2 A/ (p 2 , m 1; m 3 ) 

- 2m 2 / (0, mi, mi) - 2m 2 / (0, m 3 , m 3 ) - 2p 2 /3 



+ cos (a — P) ( [2m 2 + 2m 3 - p ]J (p , m 2 , m 3 ) - [m 2 - m 3 ] A7 (p , m 2 , m 3 ) 



2m 2 /o(0, m 2 , m 2 ) - 2m 2 / (0, m 3 , m 3 ) - 2p 2 /3 



(58) 



(P 2 ) 



a 
4ir 



^fo 2 ( 2p2 -m 2 f ,)I Q (p 2 ,mf,mfi 



2my/o(0, m/, m/) — 2m^I (0, my, my) — (m^ — m^) A/ (p , m./, m//) —2p/3 
2 



1 



12s 2 



(5p 2 + 2M^)I (p 2 , M w , A) - 2M^/ (0, M w , M w ) - M^AJ (p 2 , M w , A) + p 2 /3 
[40c 2 - l]p 2 + [16c 2 + 54 - 10/c 2 ]M^V (p 2 , M w , M z ) 



- (16c 2 + 2)M^/ (0, M w , M w ) - (16c 2 + 2)Mf/ (0, M z , M z ) 
+ (8c 2 - 2)p 2 /3 - (1 + 8c 2 )(M^ - M 2 z ) 2 AI (p 2 , M w , M z ) 

- cos 2 (a -P) ([2ml -10M^-p 2 ]I (p 2 ,M w , mi )- [ml -M^] 2 AI (p 2 ,M w , mi ) 

- 2m 2 / (0, mi, m x ) - 2M^/ (0, M w , M w ) - 2p 2 /3^j 

+ sin 2 (a - P) ( [2m 2 - 10M^ - p 2 ]/ (p 2 , M w , m 2 ) - [m 2 - M^] 2 A/ (p 2 , M w , m 2 ) 
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2m 2 2 I (0, m 2 , m 2 ) - 2M^/ (0, M w , M w ) - 2p 2 /3 



+ sin (a — (3) \[2m 1 + 2m + — p ]Iq(jp , m + , mi) — [m[ — m + ] AI (p , m + , mi) 
— 2m 2 / (0, mi, mi) — 2m+/ (0, m + , m + ) — 2p 2 /3^ 

+ cos 2 (a - 0) ( [2m 2 2 + 2m 2 - p 2 ]/ (p 2 , m+, m 2 ) - [m| - m 2 } 2 AI (p 2 , m+, m 2 ) 



— 2m 2 / (0, m 2 , m 2 ) — 2m+io(0, m + , m + ) — 2p 2 /3 
+ {[2ml + 2m + ~ V 2 ]I(s{p 2 i m + , m 3 ) - [m 2 - m 2 + ] 2 AI (p 2 , m+, m 3 ) 

- 2m3/ (0,m3,m 3 ) - 2m+/ (0,m + ,m + ) - 2p 2 /3^) |. 

The sum £) f is over all fermions in the theory, and /' is the isospin partner of /. We 
have neglected quark mixing in the W self energy. 

In the calculation of the Z width, we need in addition the various vertex correc- 
tions. Here we tabulate only those corrections involving exchange of physical scalars. 
For the full calculation, one must also include diagrams involving W and Z exchange, 
which are present also in the minimal standard model. For the pure standard model 



corrections, we refer the reader to Refs. [13, 15] . Here we give the renormalized vertex 
corrections to Z — ► bb from physical Higgs exchange. Where no ambiguity is involved, 
the arguments of groups of tensor integrals that depend on the same variables have 
been specified only once at the end of the group, 
(a) Charged Higgs vertex correction: 



+ 



a f m 2 cot 2 (3 ( 1 — 75 \ / s 2 — c 



sc 



-K 00 (M z ,m t ,m + ,m H 



TYl 

+ [v t - a t ][2K m - 1/2 + M 2 Z K U - m 2 t K ]{M 2 z ,m + ,m t ,m t ) - ^-K (M 2 z ,m + ,m u m t ) 



[v b + a b ]Ii(0,m t} m^ 
m 2 tan 2 (3 ( 1 + 75 \ / s 2 - c 2 



M? 



w 



sc 



K 00 (M 2 ,m t ,m + ,m + ) + [v t + a t ][2K 00 - 1/2 



+ M 2 K 12 - m 2 K }(M 2 ,m + ,m t ,m t ) + -±K (M 2 ,m t ,m + ,m + ) 

2sc 



- [v b - a b ]li(0,m t ,m p ) 
(b) Neutral Higgs vertex correction: 



(60) 



ST,. 



a m 



b \ ~* 



8ns 2 2M 2 V ^ 



l±7s 



2 sec/3 



sc 



tan (3 cos a sin(/? — a)K 00 (M z , 0, mi, m 3 ) 
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+ tan (3 sin a cos(/3 — a)K 00 (M z , 0, m 2 , m 3 ) + cos a cos(/3 — a)K 00 (M z , 0, mi, Mz) 
— sin a sin(/9 — a)K 00 (M z , 0, m2, M^)^ 

(COS^ G! 
— — -(2K o - 1/2 + M 2 Z K 12 )(M 2 Z: mi, 0, 0) 
COS 

gin ^ rv 

+ ^^(2^oo - 1/2 + M 2 z K 12 )(Ml m 2 , 0, 0) + tan 2 /3(2^ 00 - 1/2 + Mf Xi 2 ) (M|, m 3 , 
cos^p 

- (v b ^a b )(^ph(0,m 1: 0) 
+ ^a^A(0, m 2 , 0) + tan 2 Ph(0, m 3 , 0) + J^O, M z , 0)) 1 ) 
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